Full counting statistics in disordered graphene at Dirac point: 
Prom ballistics to diffusion 



O 

(N 

C 

1—5 



A. Schuessler, 1 P. M. Ostrovsky, 1 ' 2 I. V. Gornyi, 1 ' 3 ' 4 and A. D. Mirlin 1 ' 4 ' 5 ' 6 

1 Institut fur Nanotechnologie, Karlsruher Institut fur Technologie, 76131 Karlsruhe, Germany. 
L. D. Landau Institute for Theoretical Physics RAS, 119334 Moscow, Russia. 
A. F. loffe Physico-Technical Institute, 194021 St. Petersburg, Russia. 
"* DFG Center for Functional Nano structures, Karlsruher Institut fur Technologie, 16131 Karlsruhe, Germany. 
5 Inst, fur Theorie der Kondensierten Materie, Karlsruher Institut fur Technologie, 76131 Karlsruhe, Germany. 
6 Petersburg Nuclear Physics Institute, 188300 St. Petersburg, Russia. 

The full counting statistics of the charge transport through an undoped graphene sheet in the 
presence of smooth disorder is studied. At the Dirac point both in clean and diffusive limits, 
transport properties of a graphene sample are described by the universal Dorokhov distribution 
of transmission probabilities. In the crossover regime, deviations from universality occur which 
can be studied analytically both on ballistic and diffusive sides. In the ballistic regime, we use a 
diagrammatic technique with matrix Green functions. For a diffusive system, the sigma model is 
applied. Our results are in good agreement with recent numerical simulations of electron transport 
in disordered graphene. 
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I. INTRODUCTION 



Electron transport in graphene remains a field of in- 
tense experimental and theoretical activity^. The hall- 
mark of graphene is the massless Dirac character of low- 
energy electron excitations. This gives rise to remarkable 
physical properties of this system distinguishing it from 
conventional two-dimensional metals. The most remark- 
able effects arise when the chemical potential is located 
in a close vicinity of the neutrality (Dirac) point. In par- 
ticular, a short- and- wide sample (with width W much 
exceeding the length L) of clean graphene exhibits at 
the Dirac point pseudo-diffusive charge transport, 3 with 
"conductivity" 4e 2 /wh and with counting statistics (char- 
acterizing fluctuations of current) equivalent to that of a 
diffusive wire4~— In particular, the Fano factor (the shot 
noise power divided by the current) takes the universal 
value F — 1/3 that coincides with the well-known result 
for a diffusive metallic wire£ This is at odds with usual 
clean metallic systems, where the conductance (rather 
than conductivity) is independent of L and the shot noise 
is absent (F = 0). The reason behind these remarkable 
peculiarities of transport in clean graphene at the Dirac 
point is linearly vanishing density of states. This im- 
plies that the current is mediated by evanescent rather 
than propagating modes. The above theoretical predic- 
tions have been confirmed in measurements of conduc- 
tance and noise in ballistic graphene flakes^— Recent 
advances in preparation and transport studies of sus- 
pended graphene samples also indicate that the system 
may be in the ballistic regimei 12 ' 13 

Effects of impurities on transport properties of graph- 
ene are highly unusual as well. In contrast to conven- 
tional metals, ballistic graphene near the Dirac point con- 
ducts better when potential impurities are added;" - — 
Quantum interference in disordered graphene is also 
highly peculiar due to the Dirac nature of the carri- 



ers. In particular, in the absence of intervalley scatter- 
ing, the minimal conductivity^ ~ e 2 /h is "topologically 
protected" from quantum localization^ The exact value 
of the conductivity of such a system at the Dirac point 
depends on the type of intravalley scattering (random 
scalar or vector potential, or random mass, or their com- 
bination). For the case of random potential only (which 
is experimentally realized by charged scatterers) the con- 
ductivity in fact increases logarithmically with the length 
L, in view of antilocalizationi 14 i 19 ~ — 

In our previous worki£, we have studied the evolution 
of conductance of a short-and-wide graphene sample from 
the ballistic to the diffusive regime. We have also shown 
that the leading disorder-induced correction to the noise 
and full counting statistics in the ballistic regime is com- 
pletely governed by the renormalization of the conduc- 
tance. This implies, in particular, that the Fano fac- 
tor 1/3 remains unaffected to this order. Indeed, the 
experiment s 10 ^ 3 give Fano factor values in the vicinity 
of 1/3 at the Dirac point for different system lengths L. 
One could thus ask whether deviations from this value 
should be expected at all. 

In this work we present a detailed analysis of the shot 
noise and the full counting statistics in samples with 
long-range (no valley mixing) disorder. We show that 
to second order in the disorder strength a correction to 
the universal counting statistics of the ballistic graphene 
does arise. We calculate this correction and demonstrate 
that it suppresses the Fano factor below the value 1/3. 
For the case of random scalar potential, we also ana- 
lyze the opposite limit of large L when the system is 
deep in the diffusive regime. Generalizing the analysis 
of weak-localization effects on the counting statistics by 
Nazarov2i, we find that the Fano factor returns to the 
value of 1/3 from below with increasing L. The approach 
to 1 /3 is however logarithmically slow, These results com- 
pare well with recent numerical work o 20 ' 21 and particu- 



2 



larly with the most detailed study by Tworzydlo et al££ 
The structure of the paper is as follows. In the Sec. 
HI! we describe the general matrix Green function for- 
malism and its application to the problem of full count- 
ing statistics. The model for graphene setup and dis- 
order is introduced in Sec. IIII1 We proceed with ap- 
plying matrix Green function method to the calculation 
of the distribution of transmission probabilities of the 
clean graphene sample in Sec. |IV] In Sec. [V] we evalu- 
ate perturbative disorder corrections to the full counting 
statistics in ballistic regime. Diffusive transport through 
disordered graphene is considered in Sec. I VII within the 
sigma-model approach. The paper is concluded by Sec. 
IVIII summarizing the main results. Technical details of 
the calculation are presented in three appendices. 



II. MATRIX GREEN FUNCTION FORMALISM 

We begin with the general presentation of the matrix 
Green function approach to the full counting statistics 
of a quasi-one-dimensional system. This formalism was 
developed by Nazarov in Ref. [25l 

Consider a quasi-one-dimensional sample attached to 
two perfect metallic leads. Transport characteristics of 
the system are encoded in the matrix of transmission 
amplitudes t mn , where the indices enumerate conduct- 
ing channels (quantized transverse modes) in the leads. 
Eigenvalues of the matrix tH determine transmission 
probabilities of the system (we use the "hat" notation for 
matrices in the space of channels). Our main goal is to 
calculate the distribution of these transmission probabil- 
ities. The full counting statistics of the charge transport 
is given by the moments of this distribution or, equiva- 
lently, by the distribution itself. The first two moments 
of the transferred charge provide the conductance (by 
Landauer formula) and the Fano factor 



G= — Trtfi, 
h 



TrM 



(1) 



The starting point of our consideration is the relation 
between the matrix of transmission probabilities and the 
Green function of the system, 



t'tnn ^ \p^mVnG rfln {x ^ X ), 



(2) 



Here v m ,n are velocities in the mth and nth channel. The 
Green function is taken in the mixed representation with 
real-space coordinates in the longitudinal direction and 
channel indices in transverse direction. The positions x 
and x' should be taken in the left and right lead respec- 
tively in order to obtain the transmission matrix of the 
full system. The conjugate matrix is related to the 
retarded Green function by a similar identity. 

The Green functions are defined in the standard way 

as 



with energy e and Hamiltonian H, the latter being an 
operator acting both on x and in the channel space. 

With the help of Eq. ([2]) we can express all the mo- 
ments of transmission distribution in terms of the Green 
functions, 



Tr (*t*)" = Tr [vG A (x, x')vG R {x', x)] ' 



(4) 



where v is the velocity operator and x and x' lie in the 
left and right lead respectively. For the first moment, 
n = l, the above identity establishes an equivalence of 
the Landauer and Kubo representations for conductance. 

The complete statistics of the transmission eigenvalues 
can be represented by the generating function 



n=l 



(5) 



Once this function is known, all the moments of transmis- 
sion distribution are easy to obtain by expanding the gen- 
erating function in series at z = 0. An efficient method 
yielding the whole generating function was proposed in 
Ref.l25l It amounts to calculating the matrix Green func- 
tion defined by the following equation: 



e-H + iO 

-yfzv8(x - Xr) 



</zv5(x — xl) 
e-H-iO 



G(x,x') 

--5{x-x')i. (6) 



The parameter z here corresponds to the source field 
mixing retarded and advanced components of the ma- 
trix Green function. The positions Xl and xr, where the 
source field is applied, lie within the left and right lead 
respectively. We will refer to this specific matrix struc- 
ture as the RA (retarded - advanced) space and denote 
such a matrices with the "check" notation. 

The main advantage of the matrix Green function de- 
fined by Eq. ^ is the following concise expression for 
the generating function of transmission probabilities: 



\ z 




v 



G{x R ,x R ) 



:Tr 



v 




G(x L ,x L ) 



(7) 



(e - H ± i0i)G ' (x, x 1 ) = 8{x - x')t, 



(3) 



The validity of this equation can be directly checked by 
expanding the Green function in powers of z with the 
help of perturbation theory and comparing this expan- 
sion termwise with Eq. The equivalence of these two 
expansions is provided by the identity 

Another and, probably, most intuitive representation 
of the full counting statistics is given by the distribution 
function of transmission probabilities P(T). This func- 
tion takes its simplest form when expressed in terms of 
the parameter A related to the transmission probability 
by T = 1/ cosh 2 A. In terms of A the probability den- 
sity is defined by the identity P(T)dT = P(X)dX. The 
definition of the generating function, Eq. ([3]), implies a 
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trace involving all transmission probabilities. With the 
distribution function of these probabilities we can express 
J-(z) by the integral 



Hz] 



P(X)d\ 
cosh 2 A — 



(8) 



The function F{z) has a branch cut discontinuity in the 
complex z plane running from 1 to +oo. The jump of 
the function across the branch cut determines the distri- 
bution function fsee Ref. UM for derivation! : 



statistics of a disordered graphene sample. Our strat- 
egy is as follows. First, we calculate the matrix Green 
function of a clean rectangular graphene sample and ob- 
tain the full counting statistics with the help of Eq. (J7J . 
Then we introduce disorder in the model perturbatively. 
Evaluation of the free energy by diagrammatic methods 
yields disorder corrections to the full counting statistics 
of a clean sample. 



III. MODEL 



P(X) = [JYcosh A + iO) - JYcosh A - iO)l . (9) 

In other words, Eq. ((8]) solves the Riemann-Hilbert prob- 
lem defined by Eq. ((9|). 

The generating function F(z) can be related to the 
"free energy" of the system in the "external" source field. 
The free energy is defined in terms of the functional de- 
terminant 



Q = TrlnG, 



T 



dz 



(10) 



The free energy can be calculated by standard diagram- 
matic methods and hence provides a very convenient rep- 
resentation of the full counting statistics. It is also con- 
venient to parametrize the argument of the free energy 
by the angle according to z = sin 2 (0/2). 

Thus we have three equivalent representations of the 
full counting statistics by the functions F(z), f (A), and 
0(0) . In this paper we calculate the transport character- 
istics of a disordered graphene sample in terms of its free 
energy 51(0). The two other functions can be found with 
the help of identities 



Hz) = -?-**■ 

sin 90 
P(A) = 2 R eg 

TT O0 



-2 arcsin \fz 



(11) 

(12) 



The first of these relations directly follows from Eq. (fTOj) 
while the second one is the result of the substitution of 
Eq. (HJ) into Eq. ©. 

The two most experimentally relevant quantities con- 
tained in the full counting statistics, namely, conductance 
and Fano factor, Eq. ([1]), can be expressed in terms of 
any of the functions introduced above. Then the follow- 
ing expressions for the conductance and the Fano factor 
hold: 



_ 2e 2 d 2 n 

h 90 2 



6=0 



1 _ 2 9 4 O/90 4 
~ 3 ~ 3 d 2 n/d(f> 2 



(13) 
(14) 



6=0 



We will apply the matrix Green function formalism 
outlined in this section to the problem of full counting 



We will adopt the single-valley model of graphene. 
More specifically, we will consider scattering of elec- 
trons only within a single valley and neglect interval- 
ley scattering events. Indeed, a number of experimen- 
tal results show that in many graphene samples the 
dominant disorder scatters electrons within the same 
valley. First, this disorder model is supported by the 
odd-integer quantization ^ 17 ? 26 of the Hall conductivity, 
a xy = (2n + l)2e 2 //i, representing a direct evidence 2 ^ in 
favor of smooth disorder which does not mix the val- 
leys. The analysis of weak localization also corrobo- 
rates the dominance of intra- valley scattering 2 ^. Further- 
more, the observation of the linear density dependence 1 
of graphene conductivity away from the Dirac point can 
be explained if one assumes that the relevant disorder is 
due to charged impurities and/or ripplesi 19 ' 29 " — Due to 
the long-range character of these types of disorder, the 
intervalley scattering amplitudes are strongly suppressed 
and will be neglected in our treatment. Finally, appar- 
ent absence of localization at the Dirac point down to 
very low temperature o 17 ? 26 ' 33 points to some special sym- 
metry of disorder. One realistic candidate model is the 
long-range randomness which does not scatter between 
valley a 18 ' 34 . 

The single-valley massless Dirac Hamiltonian of elec- 
trons in graphene has the form (see, e.g., Ref. 0) 

H = v ap + V(x,y), V(x, y) = a^V^x, y). (15) 

Here <r M (with fi = 0, x, y, z) are Pauli matrices acting on 
the electron pseudospin degree of freedom corresponding 
to the sublattice structure of the honeycomb lattice, a = 
{<j x ,Oy\, and the Fermi velocity is vq ps 10 s cm/s. The 
random part V{x, y) is in general a 2 x 2 matrix in the 
sublattice space. Below we set K = 1 and Vq = 1 for 
convenience. 

We will calculate transport properties of a rectangu- 
lar graphene sample with the dimensions L x W. The 
contacts are attached to the two sides of the width W 
separated by the distance L. We fix the x axis in the 
direction of current, Fig. [TJ with the contacts placed at 
x = and x = L. We assume W L, which allows us to 
neglect the boundary effects related to the edges of the 
sample that are parallel to the x axis (at y — ±W/2). 

Following Ref. 0, metallic contacts are modelled as 
highly doped graphene regions described by the same 
Hamiltonian (1151) . In other words, we assume that the 
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FIG. 1: (Color online) Schematic setup for two-terminal 
transport measurements. Graphene sample of dimensions 
L x W is placed between two parallel contacts. We assume 
W S> L throughout the paper. 



chemical potential Ep in the contacts is shifted far from 
the Dirac point. In particular, Ep 3> e, where e is the 
chemical potential inside the graphene sample counted 
from the Dirac point. (All our results are independent 
of the sign of energy, thus we assume e > throughout 
the paper.) We also assume zero temperature, that is 
justified provided TL <C 1. 

With the boundary conditions specified above, we are 
able to calculate explicitly the matrix Green function © 
for a clean graphene sample [V(x,y) = 0] at zero energy. 
This calculation is outlined in Appendix lAl [see Eq. (|A7j) ]. 
Using this Green function, we will study disorder effects 
in the framework of the diagrammatic technique for the 
free energy. 



IV. ELECTRON TRANSPORT IN CLEAN 
GRAPHENE 

In this section we apply the matrix Green function for- 
malism developed in Sec. [IT] to the case of clean graphene. 
These results will play the role of the zeroth approxima- 
tion for our perturbation theory. 

The matrix Green function is derived in Appendix [A] 
The generating function for the full counting statistics 
is given by Eq. ((TJ). With the Green function (|A7j) . we 
obtain 



>in 2 f) = 



W 



Tr 



G(0,0;0) 



W 4> 

ttL sin t 



The corresponding dependence of the free energy on the 
source field <f> follows from integration of Eq. ([TTT) . This 
yields a simple quadratic function 



AttL 



(17) 




FIG. 2: Lowest energy correction to the free energy of the 
system. 



characteristic for disordered metallic wires 



Po(T) 



W 



1 



2ttL T\f\~^T 



(18) 



Hence electron transport in clean graphene at the Dirac 
point is often called pseudodiffusive. 

Let us now calculate an energy correction to the pseu- 
dodiffusive transport regime. In the vicinity of the Dirac 
point (we assume tL <C 1), we can account for finite 
energy e by means of perturbation theory. The linear 
term is absent due to particlc-holc symmetry of the Dirac 
point. The lowest non- vanishing correction appears in 
the e 2 order and is given by the single diagram Fig. [3J 



a = 



We 2 



dxdx' I dyTrG(x,x';y)G(x',x;—y). 



(19) 

This integral of the product of two Green functions is 
calculated in Appendix [Bj The result takes the form 



n W (eL) 2 8 
7TL sin $ 9i 







2tt 



2- 



where ip is the digamma function. 

As explained above, the free energy f2 e contains infor- 
mation about the full counting statistics, i.e. all moments 
of the transfered charge. In particular, from Eqs. (fT5|) 
and ([Ft]) we obtain the following results for the conduc- 
tance and Fano factor: 



G 



C2 



^^[l + Cl ( e L) 2 ], F =I[l + c 2 ( £ L) 2 ] 



0.101, 



35C(3) 124C(5) 

37T 2 7T 4 

28C(3) 434((5) + 4572C(7) 

1571" 2 7T 4 7T 6 



(21) 
(22) 

-0.052. (23) 



These expressions coincide with the results of Ref. [l6| ob- 
tained within an alternative (transfer-matrix) approach. 



This remarkably simple result reveals the convenience of 
the source field parametrization z — sin 2 ((/>/2). The clean 
sample responds linearly to the "external" field cf>. The 
distribution of transmission probabilities given by Eq. 
(fT2")l is just a constant, Pq(X) = W/irL, in terms of A. 
This means the distribution acquires the Dorokhov form 4 



V. DISORDERED GRAPHENE: BALLISTIC 
LIMIT 

Let us now include the random part V(x, y) of the 
Hamiltonian (|15[) into consideration. There are in total 
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For the four possible disorder types, this yields 



(a) 



FIG. 3: Loop diagrams for the disorder corrections to the 
ground state energy fl, (a) first-order, (b) and (c) second- 
order. 



four different types of disorder within the single valley 
Dirac model: Vq is the random potential (charged im- 
purities in the substrate), V x and V v correspond to the 
random vector potential (e.g. long-range corrugations), 
and V z is the random mass. We will assume the standard 
Gaussian type of disorder characterized by the correla- 
tion function 



(T^(r)K(r')) = 2^,^(11- -r'|) 



(24) 



The functions (r) depend only on the relative distance 
r and are strongly peaked near r = 0. Thus we deal 
with isotropic and nearly white-noise disorder. However, 
in order to accurately treat ultraviolet divergencies aris- 
ing in our calculation, we keep a small but finite disor- 
der correlation length. The results will be expressed in 
terms of four integral constants characterizing the disor- 
der strength 



drWftQrl). 



(25) 



Within the specified Gaussian disorder model, pertur- 
bative corrections to the free energy are given by the 
loop diagrams. The first and second order corrections 
are shown in Fig. [3l Dashed lines in these diagrams de- 
note disorder correlation functions 



A. First-order correction 

The first-order correction to the free energy f2(</>) is 
given by the loop diagram containing two Green func- 
tions and one impurity line, Fig. [3ji. The Green func- 
tion at coincident points diverges. That is why we keep 
a finite correlation length calculating the first-order di- 
agram. Assuming the separation between two vertices 
of the diagram is given by the vector d, we obtain the 
expression 

fia= / dS^w^SWP, (26) 
<SnM = 7T / dr Tr L.Gfr, r + S)a„G(r + 8, r)l . (27) 



Now we substitute the Green function from Eq. (|A7[> and 
expand the correction to the free energy in powers of S. 



n(o).w 



ttW 
2L 2 



L 






f 'I A 





2sin 2 f; \ 



6S 2 ir 2 5 2 



o(*),(i/) _ *W 



L 




Jd X 


[ -l A 





4 sin 2 2^ \ 



1 , 



(28) 



(29) 



In these expressions we encounter two types of divergent 
terms: one with negative power of S and one with an 
integral of sin~ 2 (7riz;/L), which diverges at x — and 
x = L. These terms, however, are free of the source 
parameter <f> and hence do not change any observables. 
The </)-dependent terms are finite and, after integrating 
over S in Eq. (pM]). yield the simple result^ 



Q a = const + (a<o — a z ) 



AttL 



(30) 



It provides a linear (in a M ) correction to the free energy 
of the clean sample, Eq. (|17[) . merely changing the over- 
all prefactor (conductance) but preserving the quadratic 
dependence on <j) and hence the form of the Dorokhov 
distribution. Thus the linear disorder correction does 
not destroy the pseudodiffusive character of transport in 
graphene at the Dirac point. 



B. Second-order corrections 



Since the lowest disorder correction (|30|) preserves the 
form of the Dorokhov distribution, we proceed with 
higher order corrections. Our aim is to find a deviation 
from the pseudodiffusive transport. The second-order 
correction to the free energy is due to the diagrams Fig. 
[3Jd and(3Ji. The diagram with parallel impurity lines (Fig. 
|3b) yields 



Cl b = f d5dS'Y,^(\S\)M\*\)rtr\ (3 1 ) 

n^ u) = 2tt 2 J dr dr' Tr [^Gfr, r + S)a fi G(r + S, r') 
x a u G(r', r' + S')a v G{r' + 5', r)l . (32) 



Using the Green function from Eq. (|A7|) . we expand the 
correction to the free energy in powers of d and 6' . Then 
we drop all ^-independent terms and average with re- 
spect to the directions of 6 and 6'. The following four 
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contributions to the free energy are non zero: 



(°°) - c>( zz ) - _n( 0z ) - n( z °) 

"b — "b — "b — "b 
2 pL 



1 



— — r / da; dx' 

64 Wo „ v 

1 



dy 



2 7r(x-\-x'+iy) 



2L 



2 ^[x — x'+iy) 



+ C.C. 



2L 



(33) 



After integrating with respect to x and a;' the above ex- 
pression vanishes. Thus we conclude that the diagram 
Fig. [2t> gives no contribution to the free energy, 



n h = o. 



(34) 



Let us now consider the diagram Fig. [5J; with crossed 
impurity lines. This diagram contains no Green functions 
at coincident points and hence does not require regular- 
ization. We can replace disorder correlation functions 
by equivalent delta-functions and obtain 



ft V 



Q ( f u) =tt 2 J drdr'Tr [a„G(r,r')a v G(r' , 



(35) 



(36) 



With the Green function (| A7|) we find the following con- 
tribution to the sum in Eq. ([33]) : 



n (oo) = n (oz) = n w 

L 



TT W 

64L 4 



dx dx' 



1 



24>y 

dy cosh — — - 



1 



(37) 



sin 



tt(x-\-x' +iy) I 



= n (00) + 



2L 



8L 4 



sm 



7r(:r — x' +iy) I 



2L 



■DC 



tyy 

L 



. ir(x + x' + iy) . ir(x — x' + iy) 
sm sm - 



2L 



2L 



(38) 



Two-dimensional integrals with respect to x and x' are 
straightforward due to periodicity of the integrand. As a 
result, the free energy is expressed as a single y integral: 



TT 2 W 

8L 2 



dy 



cosh(20y/L) 
sinh 2 (iry / L) 



(ato + a z ) 2 coth 



(ato + 3a z )(ao - a z ) 



■ny_ 
L 



(39) 



This integral diverges at y = 0. Expanding near this 
point, we find that the integrand behaves as (L/|y|) 3 + 
2<fi 2 L/\y\. The most singular part is ^-independent and 
hence unobservable. Integral of the second term diverges 



logarithmically and multiplies <fi 2 . This gives a logarith- 
mic correction to the conductance of the system preserv- 
ing the pseudodiffusive form of the transmission distri- 
bution. Let us cut off the logarithmic integral at some 
ultraviolet scale y = a that is the smallest scale where 
the massless Dirac model with Gaussian white-noise dis- 
order applies, e.g. the scale of the disorder correlation 
length or lattice spacing in graphene. The upper cut-off 
is already embedded in the integrand of Eq. (|39|) : the 
small y expansion is valid for y < L. Thus we can isolate 
the divergent part of the integral ([3TR) and the remaining 
Q c correction, which has a non-trivial dependence on cj). 



W<p 2 
AttL 



{(a + a z ) 2 [21n( J L/a)+ Wl (^)] 

+ (ao + 3a 2 )(a -a z )w 2 (</>)}- (40) 



Since the logarithmic term in the free energy contains an 
ultraviolet parameter a defined up to a model-dependent 
constant, the functions ui^ist 1 ) are fixed up to an arbi- 
trary constant. With this accuracy, we find 



Wl(0) 



cosh^ 

dy 



sinh 3 |^| 



cosh 



2<py 



20V 



2L(j) 2 

— OO 

const — ip((j)/7r) — ip(— <£/7r), 

OO 

3 r cosh(2(j)y/L) - 1 



L 



2L(j) 2 
const + 7r 



^ sinh 2 (ny / L) 
2 (j) cot <f) — 1 



(41) 



(42) 



The logarithmic correction in Eq. (j40|) can be included 
into an effective L-dependence of the disorder strength 
parameters a M by renormalization group (RG) methods. 
The model of two-dimensional massless Dirac fermions 
subject to Gaussian disorder and its logarithmic renor- 
malization appeared in various contexts. In particular, 
disorder renormalization in graphene was considered in 
Refs. 16.32.36. One- loop RG equations for effective dis- 
order couplings as functions of a running scale A are 



dao 
d In A 

da x 



2(a + a z )(a + a x + a y ), 



da*, 



— 2a a z 



din A dlnA 

da z 

QY n \ = ( a ° + az )\ az + a x + aty). 



(43a) 
(43b) 
(43c) 



Parameters defined in Eq. (|25[) serve as initial conditions 
for the RG equations at an ultraviolet scale a. Integrating 
Eqs. (|4"3"j) up to the largest scale that is the system size 
L, we obtain effective disorder couplings — a^L) and 
automatically take into account all leading logarithmic 
contributions like the one in Eq. f|40|) . This allows us 
to replace the disorder parameters in the free energy by 
their renormalized values and drop the logarithm from 
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FIG. 4: Functions pi and P2 entering the disorder correction 
to the distribution of transmission probabilities (|45[) . In the 
case of random scalar potential (qo), the distribution is de- 
termined by the sum pi + P2 only. 



Eq. (|40l) . Collecting in this way the contributions (fTTj) . 
(f3"0")) . and PO]) . we obtain the final expression for the free 
energy up to the second order in renormalized disorder 
parameters, 



n = 



W(f> 2 
AttL 



1 + do - a z + (a + a z ) 2 wi(0) 

+ {a + 3a z )(& -& z )u 2 ((j)) . (44) 



Thus we have established a deviation from pseudodif- 
fusive transport regime (Q ~ </> 2 ) in the second order in 
disorder strength. 



Corrections to the distribution function 



The functions p\ and p 2 cannot be used for direct calcu- 
lation of transmission moments due to their divergence 
at A = 0. This divergence signifies the break down of per- 
turbative expansion in small values of disorder couplings 
close to A = (that is T = 1). Comparing disorder 
correction with the distribution in the clean sample, we 
conclude that the result (|45[) is valid provided A > a. 

The deviation from pseudodiffusive transport regime 
can be experimentally demonstrated as a correction to 
the Fano factor F — 1/3 characteristic to the diffusive 
systems. Divergence of the functions pi.2 at A = 1 pre- 
vents us from calculating transmission moments from the 
distribution function (|45[) . However, we can obtain trans- 
port characteristics from the free energy (|44|) instead. 
With the help of Eq. (fl"4|) . we find the Fano factor up to 
quadratic terms in the renormalized disorder strength, 

„ 1 16C(3),_ „ , 2 8tt 2 ._ 

F = ~ 5— (a + a z ) + —— (a + 3a z )(a - a z ) 

6 7T Z 45 

fa ~ -0.194 &l- 0.388 a &z- 7.212 d% (47) 
o 

Remarkably, any weak disorder, irrespective of its ma- 
trix structure, suppresses the Fano factor. [Note that 
the energy correction (|2Tj) is also negative.] The correc- 
tion to the Fano factor increases with increasing sample 
length L due to renormalization (|43|) . At some length I, 
referred to as the mean free path, one of the renormalized 
disorder couplings reach a value of order unity and the 
perturbative RG treatment breaks down. This signifies 
the crossover from ballistic to diffusive transport regime. 
Disorder correction to the Fano factor becomes strong in 
this crossover region. To go beyond the mean free path 
scale we resort to other methods designed for diffusive 
systems. 



Let us now derive a correction to the Dorokhov distri- 
bution function of transmission probabilities. In the A 
representation, the distribution function is given by Eq. 
([1~2"T). Using the result (|33]l. we obtain 

W r 9 
P(A) = — 1 + &q - & z + (a + &z) Pi(X) 
ttL l 

+ (So + 3a z )(a - a z )p 2 {\) ■ (45) 

Similarly to cJi.,2, the functions Pi^A) are defined up to 
a model-dependent constant. From Eqs. (14"TT) and ((32) 
we obtain 



d 

Pi (A) = const — 2 Re — — 
o\ 

^2 

Vi W = const + 



7T 



2sinh z (2A) 



(46a) 
(46b) 



The functions p\ andp2 are shown in Fig. [4] (When the 
only disorder is ao, correction to the distribution func- 
tion is given by the sum p\ +p 2 also shown in the figure.) 



VI. DISORDERED GRAPHENE: DIFFUSIVE 
LIMIT 

When the system size exceeds the mean free path, the 
sample exhibits diffusive electron transport. On a semi- 
classical level, the system can be characterized by its 
conductivity per square in this limit. At the ballistics- 
diffusion crossover the conductivity of graphene is close 
to the quantum value e 2 /h. This signifies strong inter- 
ference corrections to transport characteristics making 
semiclassical picture inadequate. These quantum effects 
lead to one of the four possible scenarios depending on 
the symmetry of disorder: 

(i) If the only disorder is random potential (ao), the 
system possesses time inversion symmetry H = 
<j2H T a2 and falls into symplectic symmetry class 
Al l 37 ' 38 . Quantum corrections to the conductiv- 
ity are positive, leading to good metallic properties 
(large dimensionless conductivity) at large scales. 

(ii) In the case of random vector potential a x ^ y , the 
only symmetry of the problem is chirality, H = 
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—a^Ha^, signifying the chiral unitary symmetry 
class AIII. Such disorder produces no corrections 
to the conductivity to all orders and can be effec- 
tively gauged out at zero energy^. From the point 
of view of its transport properties, the system re- 
mains effectively clean and ballistic at all scales. 

(hi) If the only disorder is random mass (a z ), the Hamil- 
tonian has a Bogolyubov - de Gennes symmetry 
H = o\H T u\ characteristic for the symmetry class 
D. Upon renormalization (|4U)) the disorder cou- 
pling gets smaller and the system becomes effec- 
tively clean. This means the absence of the mean 
free path scale and hence of the diffusive transport 
regime. 

(iv) In the generic case, when more than one disorder 
type is present and all symmetries are broken, the 
symmetry class is unitary (A) and transport prop- 
erties are the same as at the critical point of the 
quantum Hall transition. 

We will concentrate on the first case (random poten- 
tial) when the system eventually acquires a large param- 
eter - dimensionless conductivity - and can be quanti- 
tatively described by the proper effective field theory - 
sigma model of the symplectic symmetry class. Our con- 
sideration in this part of the paper is closely related to 
that of Ref. [M 

Derivation of the sigma model with the source fields 
z from Eq. (j6|) is sketched in Appendix [Cj The sym- 
plectic sigma model operates with the matrix field Q 
of the size AN x AN, where N is the number of repli- 
cas. Apart from replica space, matrix Q has retarded- 
advanced (RA) and particle-hole (PH) structures. The 
former is similar to the matrix Green function while the 
latter is introduced in order to account for time-reversal 
symmetry of the problem. We will denote Pauli matrices 
in RA space by A XlVlZ . Two constraints are imposed on 
Q, namely, Q 2 = 1 and Q — Q T . This yields the tar- 
get space Q E 0(AN)/0(2N) x 0(2N) characteristic for 
symplectic class systems. The sigma-model action is^ 



S[Q] = ^f drTr(VQ) 2 



(48) 



Here a is the dimensionless (in units e 2 /h) conductivity 
of the two-dimensional disordered system. The source 
field is incorporated into boundary conditions, 



Q\ 



x=0 



A 2 , Q\ x=L = A z cos(j> + A x ain<j). (49) 



The free energy of the system in the source field <fr is 
expressed through the N — > limit of the sigma-model 
partition function as 



fi = lim — ( 1 - f DQe- s ^) 



(50) 



In a good metallic sample with a S> 1 , the Q integral in 
Eq. (|50l) can be evaluated within the saddle point approx- 
imation. The action (|48[) is minimized by the following 




FIG. 5: Cooperon correction to the free energy in the diffusive 
limit. 



configuration of the field Q: 



= U- X A Z U, U = cxp 



•a ^ x 
V 2L 



(51) 



Replacing the integral in Eq. ([5T)]) with the value of the 
integrand at the saddle point, we obtain the semiclassical 
result for the full counting statistics, 



fin = lim 

JV->0 



S[Q ] Wa<t> 2 



N 



4L 



(52) 



This yields the Dorokhov distribution of transmission 
probabilities in diffusive two-dimensional system^ In or- 
der to find corrections to this result, we take into account 
fluctuations of the field Q near its saddle-point value Q - 
This is equivalent to the calculation of a Cooperon loop, 
Fig. [SI carried out in Ref. H3. 

Small fluctuations of Q near the saddle point Qq are 
parametrized by the matrix B as (we write expressions 
involving B up to the second order) 



Q = U- l A z 1 



B 2 



U. 



B = 



b 
-b T 



This parametrization of Q automatically fulfils the con- 
ditions Q 2 — 1 and Q = Q T . The sigma-model action 
expanded up to the second order in B takes the form 



S[Q] = S[Q ] 



16 



drTr\(VB) 2 -^{A x ,B} 2 . 

(54) 

Curly braces denote anticommutator. Let us separate 
B into the parts commuting and anticommuting with 
A x . These two parts do not couple to each other in the 
quadratic action (|54"1) and only the former one couples to 
the source parameter </>. Thus we can constraint the ma- 
trix B by requiring its commutativity with A^. In terms 
of b this yields b = — b T and the action becomes 



S[Q] = S[Q 



drTr 



V6Vfe T - ^ bb T 



(55) 



This quadratic form is diagonahzed in momentum rep- 
resentation. Component of momentum perpendicular to 
the leads takes quantized values ixnjL with positive in- 
teger n due to geometrical restrictions [boundary condi- 
tions (HHJ) fix b = at the interfaces with metallic leads] . 
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Momentum parallel to the leads is continuous and un- 
restricted. For each value of the momentum there are 
N(2N — 1) independent matrix elements in b. Calculat- 
ing the Gaussian integral in Eq. (|5H)) we obtain the free 
energy 

n = no-yf; / ^-ln(n 2 n 2 + q 2 L 2 -^ 2 ) 



n=l 



W 

2L 



yj ir 2 n 2 — cf) 2 



n=l 



(56) 



In the result (j5"6")) . the sum diverges at large n. The 
situation is similar to what we have encountered in the 
ballistic regime. Expanding the sum in powers of ef>, we 
see that the most divergent term is ^-independent, while 
the next term multiplies tfi 2 and diverges logarithmically. 
This is nothing but the weak antilocalization correction. 
It renormalizes the conductivity but does not deform the 
full counting statistics. Logarithmically divergent sum 
is cut at n ~ L/l where I is the mean free path. At 
larger values of n the diffusive approximation (gradient 
expansion in the sigma model) breaks down. In terms of 
rcnormalized conductivity, the free energy reads 



2L 



~ , 2 00 / 



1 L 1 , L 
a = a H In — m — In — . 

7T I 7T I 



2im J 



(57) 
(58) 



The bare value of conductivity, er, is of order one and 
hence negligible in comparison with the large renormal- 
izing logarithm. The sum over n in Eq. (|57p is convergent 
and provides the deviation from semiclassical Dorokhov 
statistics of transmission probabilities. 

In fact, a more rigorous procedure is to perform first a 
renormalization of the sigma model from the mean free 
path scale I to the scale ~ L. Then the free energy can 
be calculated perturbatively. It turns out, however, that 
this yields a result identical to the one obtained above 
within the perturbative analysis at the scale I. Indeed, 
the RG equation da/d In A = l/ir will lead exactly to 
the renormalization of conductivity a i-> a, see Eq. (|58[). 
The consequent evaluation of the perturbative contribu- 
tion to 51 yields Eq. ([56]) with a replaced by a and the 
sum restricted to a finite (independent of L) number of 
terms. In other words, the renormalization shifts the log- 
arithmical contribution to a from the second to the first 
term in square brackets in Eq. (f56|) . 

Let us derive the distribution function P{\) from the 
free energy (I5T)) . Applying Eq. (TT2"]) . we obtain the result 
in the form 



p(X) = ^-[a + p(X)}., 



7T 



Re- 



2i\ 



1 



V^ 2 



(n + 2i\) 2 n 



(59) 
(60) 



0.2 

0.0 : 

2 -0.2 : 
ft, 

-0.4 : 
-0.6 : 







2 3 
A 



FIG. 6: Correction to the distribution of transmission eigen- 
values in the diffusive limit. 



At small values of A, the sum in Eq. (|60p is determined 
by the term with n = 1. In the opposite limit, the sum 
can be estimated by the corresponding integral with the 
help of Euler-Maclaurin formula. Thus we obtain the 
asymptotic expressions 



P(A) 



1 

8^A' 



A< 1, 



(61) 



-- In A, A > 1. 

7T 



The function p(X) is shown in Fig. |6l It is qualitatively 
similar to the numerical result of Ref. 

Deviation from the semiclassical transport can be 
demonstrated by the correction to the Fano factor. With 
the help of Eq. ([T4]) . we obtain 



2C(3) 



I 



r3^ 



0.244 
ln(£/Z)' 



(62) 



A similar correction to the Fano factor was found numer- 
ically in Ref. [22]. We compare the numerical results with 
Eq. ([12]) below. 

In the case of weak scalar disorder (described by 
the coupling ao), the system undergoes a continuous 
crossover from ballistic to diffusive transport regime as 
the size L grows. In both limiting cases, we encounter 
nearly Dorokhov distribution of transmission probabili- 
ties with small corrections, Eqs. (jT7|) and (|62p . on both 
sides of the crossover. In the ballistic limit, we can 
formally introduce a dimensionless conductivity as a — 
(L/W)G/(e 2 /h). Then the corrections to the Fano factor 
are expressed in terms of the conductivity 



16C(3) 8 



1 

F =3 



45 ' V 



2C(3) 



l) 2 , 7R7-KI, 



a > 1. 



This Fano factor as a function of conductivity is shown 
in Fig. [7] together with numerical results from Ref. [U In 
the numerical simulations, a single valley of graphene was 
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FIG. 7: Fano factor as a function of conductivity. Solid lines 
show ballistic and diffusive results (|63l) . Dashed line corre- 
sponds to the asymptotic value F = 1/3. Solid symbols are 
numerical results from Ref. [2^, the size of rectangles corre- 
sponds to the error estimate. 

modeled using a finite-difference approach. By construc- 
tion, disorder in Ref. [22| has the symmetry of scalar po- 
tential which does not mix the valleys. It is this symme- 
try (class All) which is considered in the present section. 
Our results perfectly agree with the numerics in the dif- 
fusive limit (see Fig. [7]) in the range ira > 3. On the bal- 
listic side, the deviation is due to the non-universality of 
the ballistic transport. Specifically, the function F(a) de- 
pends crucially on the microscopic details of disorder. In 
the numerical analysis of Ref. [22], the model with strong 
scatterers was used, while in the present paper we adopt 
the model of weak Gaussian white-noise disorder. For 
theoretical predictions on electron transport in the pres- 
ence of strong scatterers see Ref. H3- 

An earlier numerical study of Ref. [2(| based on the 
transfer-matrix description of the Dirac problem, re- 
ported the value of the Fano factor in the range 0.29—0.30 
(for different samples) with the conductivity, 7tct, of the 
largest systems varying from 6 to 10. This is consistent 
with our predictions for the diffusive transport regime 



(see Fig. [7]). The behavior of F in the ballistic regime is 
different due to the reasons described above (strong vs. 
weak disorder). A non-monotonous dependence F{a) at 
the Dirac point was also observed in Ref. [2l|. 

The Fano factor is 1/3 both in the clean and strongly 
disordered limits. In the crossover from ballistics to dif- 
fusion, the Fano factor strongly deviates from this uni- 
versal value signifying the breakdown of the (pseudo-) 
diffusive description characterized by Dorokhov distribu- 
tion of transmission probabilities. 



VII. SUMMARY 

We have studied the full counting statistics of the 
charge transport through an undoped graphene sheet in 
the presence of weak and smooth (not mixing valleys) dis- 
order. We have identified deviations from the Dorokhov 
distribution of transmission probabilities both in ballistic 
[Eq. (g5j), 05])] and diffusive [Eqs. ([55 ]) , ([50 ]) ] regimes. In 
the former case, corrections are model-dependent while 
in the latter case only the symmetry of disorder matters. 
We have considered Gaussian white-noise disorder in the 
ballistic regime and potential disorder (symplectic sym- 
metry class) in diffusive limit. Deviation from (pseudo-) 
diffusive transport always results in a negative correction 
to the Fano factor, F < 1/3. Our results are in good 
agreement with recent numerical simulations of electron 
transport in disordered graphene, see Fig. [7J 
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Appendix A: Matrix Green function 



The full counting statistics of the electron transport is conveniently expressed in terms of the matrix Green function 25 
in the external counting field z = sin 2 (^/2), Eq. ([6]). For the clean graphene sample attached to perfect metallic leads, 
Fig. [TJ this Green function satisfies the following equation: 

Since the operator in the left-hand side of the above equation commutes with the y component of the momentum, 
we will first calculate the Green function in the mixed coordinate-momentum representation, G p {x,x'). Inside the 
sample this function satisfies 

G p (x,x') = S(x-x'). (A2) 



d 



11 



We will look for a general solution of this equation in the form 

G p (x,x') = e T.p(*-i/2) Afe «r.p(*'-£/a) ) 



M 



M < , x < x', 
M > , x > x'. 



(A3) 



The chemical potential profile together with the infinitesimal terms ±i0 in Eq. (|Aip defines the boundary conditions 
for the Green fi 
we thus obtain 



for the Green function. The counting field z can also be incorporated into the boundary conditions. In terms of M< 



1 1 i\fz iy/z 
1 -1 



-<t zP L/2 m _ 



0. 



1 -10 

—i\fz —i\fz 1 1 



M-> = 0. 



(A4) 



Delta function in the right-hand side of Eq. (JA2J) yields a jump of the Green function at x = x' which provides the 
relation 



Af> - M< = ~ia x . 

The matrices M<, and hence the Green function, are completely determined by Eqs. (|A4[) and (|A5 



(A5) 



2(cosh pL — z) 



z + 



coshpL z 
sinh 2pL 



2 



sinh 2pL 
2 

coshpL 



i\/ze 



-pL 



v ^ 



iwze 



-pL 



cosh pi — z — 
sinh 2pL 



*V2 \ 

sinh 2pL 



± 



(A6) 



coshpL J 



Fourier transform in p yields the Green function in the full coordinate representation. To facilitate further calcula- 
tions, we decompose this Green function into the following product of matrices: 



I i COsh sinti \ ( 



G {x,x';y) = — U(x)A 



2L Sinh 2L 



sinh 



2L 



-i cosh 



4>y 



sin yj^(x + x' + iy) sin (x — x' + iy) 
1 1 



AU^ix'), (A7) 



A = 



a z 
1 



2L' ra \sin^(x-x'-iy) sin ^ (x + x> - iy) j g 
U{x) 



4>(L — x) cj)(L — x) 1 
sin — — cos ■ 



RA 



2L 

4>x 



2L 

<j)X 

lSm 2L 



(A8) 



2L 2L ' RA 

Here we have used the source angle (f> defined by z = sin 2 (</>/2). The matrices U{x) and U~ 1 (x r ) operate in the 
retarded-advanced space only and hence commute with any disorder operators placed between the Green functions. 
As a result, factors U and U~ 1 drop from expressions for any closed diagrams. The matrices A in the above equation 
allow us to decompose the Green function into a direct product of the two operators acting in the RA space and in 
the sublattice space. 



Appendix B: Energy correction to the full counting statistics 

In this appendix we evaluate the diagram Fig. [2] for the lowest energy correction to Sl(</>). Substituting Green 
function (|A7I) into Eq. (|19[) and performing rescaling of integration variables we obtain 



a 



WLe 2 



oc 

dx dx' J dy cosh(</>y) 



1 



1 



cosh(7r?/) — cos7r(x + x 1 ) cosh(7ry) — cos7r(x — x 1 ) 



(Bl) 



The first (second) term in square brackets depends only on sum (difference) of x and x' . This allows us to integrate 
over the difference (sum) of these variables. After some shifts of variables the remaining integral takes the form 



a = - 



WLe 2 



du u sin ■ 



dy cosh(</>y) 



WLe 2 



cosh 2 (Try) - sin 2 (nu/2) 2sin(</>/2) 



du 'i 



sin(0u/2) 
cos(7ra/2) 



(B2) 
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The last expression is the result of y integration. It can be performed, e.g., by closing the integration contour and 
summing up residues in the upper half-plane of imaginary y. In order to make this sum convergent, one has to add a 
weak damping factor by an infinitesimal imaginary shift of <f>. 

We proceed with the last integral in Eq. (|B2I) by representing 1/ cos(7tm/2) as a Fourier series 



a 



2WLe 2 d 
sin(>/2) 90 



f eh °° 

I du cos^^(-l)"cos[7ru(n + l/2)]. 



(B3) 



Convergence of this Fourier series should also be justified by a proper damping factor. This does not change the final 
result of the calculation hence we omit such extra factors for simplicity. Performing the integration over u we obtain 



_ 2WLe 2 d 
e ~ sin(0/2) del C ° S 2 



IE 

n=0 



1 



1 



7r(2n+l) + 7r(2n + l) 



(B4) 



The sum over n diverges logarithmically. However, this divergence is independent of <f> and hence does not influence 
any observable quantities, which are expressed as derivatives of the free energy. We can easily get rid of the divergent 
part by subtracting a similar sum over n with = 0. This yields the final result 



2I^Le 2 d 
e ~ sm(0/2)90 C ° S 



2 ^ 

n=0 



1 



1 



7r(2n + l) + ir(2n +!)-(/) 7r(2n+l) 



W (eL) 2 d 



ttL 



sin ■ 



1> 



2tt 



7T — 
27T 



41n2 + 27 



(B5) 



Here ^ is the digamma function and 7 is the Euler-Mascheroni constant. The last expression yields Eq. ([20]) of the 
main text (where we drop the unobservable constant). 



Appendix C: Derivation of the sigma model 



In order to carry out a parametrically controlled derivation of the sigma model, it is convenient to consider a 
modified problem with n>l flavours of Dirac fermions. To perform the disorder average of the free energy, we also 
introduce N replicas. (Alternatively, one can use supersymmetry. As we will treat the sigma model perturbatively, 
the two approaches are fully equivalent.) 

The derivation of the sigma model starts with the fermionic action generating the matrix Green function ([6]). 

S[<t>, 0*] = [ dr 4>? ({*0A, <rp ^a x [A+S(x) + A-6(x - L)] }s ab - V ab (r)\ ffi. (CI) 

J a,b,a ^ ' 

Here A± = (A x ± iA y )/2 are matrices operating in RA space. This action is the functional of two independent 
Grassmann two-component (in a space) vector fields <j) and ef>* . Lower indices, a and b, refer to flavours while the 
upper index a enumerates replicas. Overall, there are 4nN independent Grassmann variables in the Lagrangian. 
The random matrix V a b is symmetric, that insures the time-reversal symmetry of the model. We assume Gaussian 
white-noise statistics for the matrix V defined by the correlator 

(V ab (r)V cd (r')) = ^[6 ac S bd + S ad 6 bc ]d(r - r>). (C2) 
n 

Using the time-reversal symmetry, we rewrite the action in terms of the single four-component field tp (and its 
charge-conjugate version tp, that is linearly related to ip): 

iJ = 71 (<MV ' ^ = lJpT(7yTx = 7! z0T(Jy ) ' (C3) 

This introduces an additional particle-hole (PH) structure of the fields. Pauli matrices operating in PH space are 
denoted by t XiV _ z . Bar denotes the charge conjugation operation which has two important properties: ipi^ = ^2^1 
and (V'iV'2) T = T x o y '\bi^>\Qy't~x- The action takes the following form in terms of ip: 

S[^j] = [drJ2^a ({*0A 2 (Tp yfar x [p+S(x) + p^5(x - L)] }s ab - V ab (r)\ ^. (C4) 

3 a,b,a ^ ' 
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In this expression we have introduced the notation p± = {A x t z ± iA y )/2. 

Now we are ready to average e~ s over the Gaussian disorder distribution with the correlator (|CJ2|) . This yields an 
effective action with the quartic term. Using the above-mentioned properties of charge conjugation, we recast the 
action in the form 



dr 



a,b,a,0 



(C5) 



Next, we decouple the quartic term introducing an auxiliary 8iV x 8-ZV matrix R by the Hubbard-Stratonovich trans- 
formation. This yields the action 



S[R, i>] 



dr 



87ran 



Tri? 2 



a,a,P ^ ' 



(C6) 



Parameter 7 is an arbitrary number at this stage, its value will be fixed later. Matrix R a p couples to the product 
SaV'aV'a- This allows us to impose the corresponding symmetry constraint on the matrix R: R = a y T x R T <j v t x . 
Finally, we integrate out the fermionic fields and obtain the action operating with the matrix R only, 



S[R] = 77 Tr 2JL 
1 1 2 V47ra 



In 



iijR - erp - y/za x [p+S(x) + p-5(x - L)] X\ 



(C7) 



The bold 'Tr' symbol implies the full operator trace including integration over space coordinates. 

Derivation of the sigma model proceeds with the saddle-point analysis of the action (|CJ7[) in the absence of the 
source field z. We first look for a diagonal and spatially constant matrix R minimizing the action. The saddle-point 
equation is identical to the self-consistent Born approximation (SCBA) equation for the self energy —ijR, 



ijR = 27rao 



— ^{ijR-trp) . 



(27T) 



(C8) 



We fix 7 to be the imaginary part of the SCBA self energy, 7 = Ae _1 / Q ° with A being ultraviolet energy cut-off (band 
width). Then the saddle-point configuration for the matrix R is simply R = A 2 . This fixes the boundary conditions 
for the matrix R at the contacts. Since the leads are very good metals and fluctuations of R are strongly suppressed 
there, R = A z for x < and x > L. 

The matrix R — A z is not the only saddle point of the action (|C7|) . Other configurations minimizing the action 
can be obtained by rotations R = T~ 1 A Z T with any matrix T which commutes with erp and preserves the constraint 
R = a y T x R T cr y T x . Matrix T, and hence R, is trivial in a space. This allows us to reduce the dimension of R to AN x AN 
operating in A, r, and replicas only. The saddle manifold generated by matrices T is 0(AN)/0(2N) x 0(2N). 

Let us now restore the source term in the action and establish boundary conditions for R. The matrix R has a 
jump at the interfaces with the leads due to the delta functions in the action (|C7I) . However, we can eliminate these 
jumps by a proper gauge transformation. Let us perform a rotation R = ARA^ 1 with an ^-dependent matrix A. The 
action acquires the following form in terms of R: 



S[R] = -Tr 



~f 2 R 2 



In 



ijR - erp + io x A- x I — + i^/z[p+8{x) + pS(x - L))A \ 



(C9) 



The source field drops from this action if we choose A such that the expression in curly braces vanishes. This yields 



1, x < 0, 

A = { 1 — iy/zp + , < x < L 1 

(l - iy/zp-) (l - iy^zp + ), x > L. 



(CIO) 



Note that the matrix R, defined with the help of the above matrix A, fulfils the condition R = t x R t t x . Since delta 
functions disappear from the action, we can infer that R is continuous at the interfaces with the leads. In the left 
lead we have R = R = A z . This is the left boundary condition for the matrix R. The right boundary condition is 
fixed by the identities R = A~ X RA and R — A z for x > L. This yields 



R(L) = (1 - 2z)A z + iz 3/2 A x + 41(2 - z)A v t z . 



(Cll) 
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We can further simplify this bulky expression by performing a constant rotation R = B 1 QB with the matrix 



B 



2V2 



(1 - Z)- 1 /4 (1 + AzTz) _ i(1 _ z) l/4 (1 _ ^ 



After such a rotation the action and boundary conditions become 

~i 2 Q 2 



S[Q] = -Tr 
Q(0) = A„ 



m^Q-erp) 

Zirao 

Q(L) = A z cos cj) + A x sin ( 



(C12) 

(C13) 
(C14) 



Thus we have reduced the boundary conditions to the form (|4U1) . The matrix B is chosen such that B T B = t x . Hence 
the matrix Q obeys the symmetry constraint Q = Q T . 

The last step of the sigma-model derivation is the gradient expansion in Eq. (IC13I) . This expansion is straightforward 
for the real part of the action^ 



ReS[Q] 



■ Tr In (ijQ - erp) ( - ijQ - <rp) = - - Tr In (p 2 + j 2 + 70-VQ) 



16tt 



Tr(VQ) 2 



(C15) 



The Drude conductivity of the two-dimensional sample with n flavours of massless Dirac fermions at the Dirac point is 
(n/7r)(e 2 //i). With the dimensionless conductivity a — n/ir, we finally obtain Eq. (|48[) supplemented by the boundary 
conditions (|49p. 
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